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Abstract: We investigate neutralino dark matter in the framework of NMSSM per- 
forming a scan over its parameter space and calculating neutralino capture and an- 
nihilation rates in the Sun. We discuss the prospects of searches for neutralino dark 
matter in neutrino experiments depending on neutralino content and its main anni- 
hilation channel. We recalculate the upper limits on neutralino-proton elastic cross 
sections directly from neutrino telescopes upper bounds on annihilation rates in the 
Sun. This procedure has advantages as compared with corresponding recalcalations 
from the limits on muon flux, namely, it is independent on details of the experiment 
and the recalculation coefficients are universal for any kind of WIMP dark matter mod- 
els. We derive 90% c.l. upper limits on neutralino-proton cross sections from the results 
of the Baksan Underground Scintillator Telescope. 
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1. Introduction 

Totality of phenomena observed in neighbouring and far cosmos indicates an existence 
of unusual forms of matter and energy, which constitute up to 95% of the energy density 
of the Universe. Nowadays experiments aimed at measurements of relic anisotropy ||T|, 
gravitation lensing || and dynamics of galaxies and clasters || allow to determine 
quantitatively dark matter density with a high precision and thereby greatly developed 
early F.Zwicky hidden mass hypothesis in the Universe Q. 

In more probable cosmological model ACDM (lambda cold dark matter) the contri- 
bution of baryonic matter in total density is less than 5%, while the amount of dark mat- 
ter is about 5 times larger. Experimentally the abundance of dark matter Q x expressed 
in units of critical density is fixed in a very narrow range 0.1109 < Q x h 2 < 0.1177 [Q, 
where h is the Hubble parameter in 100 (km/s)/Mpc units. 

Basic question about nature of dark matter is still open. There are no candidates in 
the Standard Model (SM) of particle physics which could play a role of the collisionless 
dark matter in the ACDM cosmology paradigm. Among the most popular extensions of 
the SM are supersymmetry models (SUSY) which rather naturally provide with a stable 
WIMP (weakly interacting massive particle ||) — the lightest neutralino. If thermally 
produced in the early Universe the lightest neutralino undergoes self-annihilations and 
freezes out at some relic density Q x (for review see, for example, 0, §). 
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According to the electroweak (EW) theory weak interactions are mediated by gauge 
bosons while in SUSY models new particles with mass of order of EW scale also appear 
in the intermediate states. Due to the large mass of the gauge bosons weak interac- 
tions have very small interaction radius and the magnitude of the cross section of weak 
processes is a few orders smaller than that of processes mediated by strong or elec- 
tromagnetic forces. Correspondingly, annihilation cross section of WIMP dark matter 
particles ensures Q x density at EW scale. In this scenario WIMP masses are assumed 
to be between a few GeV and hundreds of TeV. 

In light of claimed spectral peculiarities the latest results of dark matter searches 
both in direct detections (DAMA [[|, CDMS II flIU| ) and in cosmic ray experiments of 
gamma, electron-positron and proton-antiproton measurements (WMAP fil l, EGRET 
H, FERMI pi, INTEGRAL 0, PAMELA 0, ATIC |6|) have been interpreted 



by many authors as an evidence of the existence of light dark matter. In particu- 
lar, neutralino dark matter in supersymmetric models were recently discussed in this 
context in papers [|17], p~8f| . 

In this note we pay attention to neutralinos in a mass range about 20 — 2000 GeV/c 2 
in the framework of Next to Minimal Supersymmetric Standard Model (NMSSM). 
In particular we consider two samples of its parameter space: mSUGRA-motivated 
models and general models without unification of parameters. With application of 
NMSSMTools package [|19|, |2(| ^TJ we perform a scan over parameter space of the 
model and calculate capture and annihilation rates for neutralinos which could be 
trapped by the Sun. We compare the prospects for NMSSM neutralino searches at 
neutrino telescopes with that of in direct detection experiments. We also discuss the 
NMSSM benchmark points [2^| in context of direct searches and data from neutrino 
telescopes. Then we consider a procedure to derive the limits on neutralino-nucleon 
cross section from the limits on annihilation rate in spirit of work |23| and use the data 
of the Baksan Underground Scintillator Telescope to derive the new limits. 



2. Theoretical frames of the model 

In supersymmetric models the neutralino is a Majorana fermion. The interactions 
of such a WIMP with ordinary matter can be both spin-dependent (SD) and spin- 
independent (SI). Moreover the dynamics of the annihilation and scattering processes 
depends on SUSY composition of the lightest neutralino, e.g. either it is a purely 
gaugino, higgsino or their mixture. We consider Z3 symmetric NMSSM model and use 



NMSSMTools package [19, 20, 21] to explore its parameter space for further analysis. 
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The relevant part of superpotential W has the form 

W = XSH U -H d + ^KS 3 , (2.1) 

where S is a superfield which is a singlet with respect to SM gauge group and H u and 
H d are two Higgs doublets, H u = (h^, H®\ , H d = (H d , H d ) ■ Here we denote 
superfields by letters with hat and their scalar components by letters without it. 
Soft SUSY breaking terms for this model read 

-C soft = m 2 Hu HlH u + m 2 Hd H\H d + m 2 s S*S (2.2) 

+ ^\A\H U ■ H d S + ^kA k S 3 + h.c.^ + l - [MtBB + M 2 W i W i + M 3 G a G a 

Here rriH u , ^H d and mj are soft scalar masses, A\, A K are trilinear coupling constants 
and M lj2 3 are soft gaugino masses. The Higgs sector of the model contains the following 
independent parameters: A, k, A K , A\, tan/3 = (H u ) / (Hd) and /i = X(S). 

In the NMSSM the lightest mass eigenstate of neutralino can be written as a linear 
combination 



X 



Z 0l B + Z 02 W 3 + Z m H° u + Zo4#2 + Z 05 s (2.3) 



of bino B, wino W 3 , two higgsinos H® and and also singlino s which is fermion 
component of the singlet superfield. Here coefficients Z^ can be chosen as real-valued. 

Below we will discuss the prospects of dark matter searches for two classes of the 
NMSSM model to be referred as Set 1 and Set 2. 

• Set 1. Here we take mSUGRA-like NMSSM with unification of soft SUSY- 
breaking parameters which are specified at GUT scale 

M lj2 , 3 = Mi/a, M~ qi = M Hi = m , M = A , (2.4) 

where are soft squarks masses and Ai are all soft trilinear coupling constants 
except for A K . We take these parameters in the following ranges: 
GeV < m < 5000 GeV, GeV < M 1/2 < 4000 GeV, 

-4000 GeV < Aq, A k < 4000 GeV, sign n > 0, 0.0 < A < 0.7, 1.7 < tan/3 < 54. 



The dark matter in the mSUGRA-like NMSSM was already discussed in Refs. [24 



25 1 . In this case A K is considered as an independent parameter and k is obtained 



from minimization of the scalar potential (see Ref. |24] for details). 

Set 2. Here we assume general NMSSM without unification of the soft terms 
and adopt the following ranges for parameters specified at SUSY breaking scale: 
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GeV < n, Mi, M 2 < 2000 GeV, 

300 GeV < M 3 < 3000 GeV, -2000 GeV < A A , A K < 2000 GeV, 
0.0 < A < 0.7, 0.0 < k < 0.6, 1.7 < tan/3 < 54; soft sfermion masses and 
sfermion trilinear terms are varied within ±3 TeV. The ranges of defined param- 
eters allow to construct rather light neutralinos. In a view of the prospects for 
the indirect dark matter searches there was related discussion in Ref. [EBf for 



m x < 100 GeV/c 2 , while we do not restrict ourselves by only light neutralinos. 
We take m t = 171.4 GeV for the mass of top-quark. For both model classes we 



scan over their parameter space by use NMSSMTools [19], ^(| pifl . All unphysical or 
phenomenologically unacceptable models are rejected by checks for absence of Landau 
pole up to GUT scale for A and k, as well for absence of unphysical global minimum 
of the scalar potential in Higgs sector and also for absence of colour breaking minima. 
Finally, we follow the NMSSMTools implications of the experimental constraints (full 



list of them can be found in Refs. [Gjl 20, Ell). Note, however, that we do not 



impose the condition that the SUSY contribution to the anomalous magnetic moment 
of muon should explain the present 3a difference between SM prediction and BNL 
results. For our further analysis we take the lightest neutralino as potential dark 
matter candidates if their relic abundance is Qdm < 0.3 and finally, we obtain about 
5 • 10 5 phenomenologically accepted models for each set of parameters presented above. 

3. Neutralino-nucleon SI and SD cross sections 

Among relevant physical quantities spin-dependent and spin-independent neutralino- 
nucleon scattering cross sections are more probable to be tested in experimental searches 
for dark matter interactions. The expected signal rate has direct dependence on neu- 
tralino scattering cross sections off ordinary nucleons not only in direct detection exper- 
iments but as well as in neutrino telescopes. In the case of annihilations of neutralino 
pairs inside the Sun where solar chemical composition contains more than 73% of hy- 



drogen |[27|| ) either SD or SI parts of the cross section could give dominant contribution 
to neutralino capture and annihilation rates. 

The effective interaction lagrangian responsible for the spin-dependent part of the 
neutralino-nuclei cross section is 

C A = d g (xY7 5 X ) {qirfq) (3.1) 

and the expression for SD cross section has the form 

32 

a SD = -G 2 F m 2 r A 2 J(J+l), (3.2) 

7T 
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where m r is the reduced neutralino mass m r = rnx . m , m and J are the mass and the 
spin of the nucleus respectively and 



A = j(a p (S p )+a n (S n )) (3.3) 



where (S p ) and (S n ) are the spin content of the proton and neutron while 

da » X ^ dn 



a = y_^L_ A W a„ = y-^A("). (3.4) 

As default values for the deltas we use 

= 0.842, A<f } = -0.427, A< p) = -0.085. (3.5) 

The axial-vector coefficients for neutron can be obtained from the values (ft.5|) by cor- 
responding isospin rotation. 

The spin-independent interactions of neutralino with nucleus are generated by the 
following effective interaction lagrangian 

£ = C Q (xx) (??) • (3-6) 

We do not discuss here the dependence of physical quantities on parameters of the 
NMSSM lagrangian and so we do not present explicit expressions for d q and c q . The 
details of the calculation of these coupling constants from NMSSM lagrangian can be 
found e.g. in Ref. [p8~| . 

The SI cross section of neutralino on nucleus can be written as 

A IT) ^ 

a sl = ^(Zf p + (A-Z)f n )\ (3.7) 

7T 

where Z and A are the charge and the atomic number of the nuclei, respectively, and 
the nucleon formfactors have the form 

— = E + ht® E wh ™ #3 = i - E < 3 - 8 > 

q=u,d,s q=c,b,t H q=u,d,s 

The parameters ftp are defined from expressions 

%/? = (N\m q qq\N) = m q B^ q N \ (3.9) 



As discussed in several recent papers |34], |35[ the main uncertainty in determination 
of SI cross section comes from the pion-nucleon scattering sigma term, cr n N, which 
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determines the coefficients fx q as follows 



h 



N 



Ba 



2a 7 



N 



'"d 



55 ) 



2B X 



y 



(3.10) 



(3.11) 



We use m u /m d = 0.553 and m s /m d = 18.9. Following Ref. we define the quantity 



-6^ — -B.o 



1.49 



and strange quark density 



where 



m u + m d 



(5„ + B d — 2B S 



and ovAr-term can be written as 



m u + m d 

&ttN = ^ \E>u + £>d) 



(3.12) 



(3.13) 



(3.14) 



(3.15) 



The recent estimates show that the value of a^N lies in the interval 55 MeV < o^n < 
73 MeV and also it is suggested that a = 35 ± 5 MeV . We use a wN = 64 MeV as 
default value for our calculations. We will back to the discussion of influence of known 
uncertainties in the nuclear matrix elements on the neutralino annihilation rates in the 
next sections. 



4. Neutralino capture rate in the Sun 

At present the estimated local density of relic dark matter in the solar system p l ° c is 
about 0.3 GeV/cm 3 . The speeds of halo WIMPs crossing the ecliptic plane are non- 
relativistic with mean velocity v ~ 270 km/s (the root-mean-square of the dispersion) 
in the Maxwell-Boltzmann distribution. As far as they are gravitationally trapped 
by the Sun, a multiple scattering of WIMP off solar matter is likely happened with 
energy losses. Due to these interactions WIMPs can be captured inside the Sun. Their 
orbital motion changes to smaller radius and toward the dense centre, where they 
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finally settle and accumulate and as well self-annihilate at a distance smaller than 
annihilation length. During the lifetime of the Sun these two processes of capture and 
annihilation may reach approximate equilibrium and in this case accumulated number 
of WIMP does not change in time. One of the prominent signature of the neutralino 
annihilations is high energy neutrinos generated in decays of annihilation products, 
since these neutrino events can be detected by neutrino telescopes. 

Such a scheme has been considered by A.Gould |30f and, e.g., later in detailed 
paper [31], see also references therein. Follow them we can obtain neutralino annihila- 
tion rate from the capture rate in the Sun by solving known evolution equation for the 
number of WIMPs 

d -^ = C c - C A N\ (4.1) 
at 

where Cc is the capture rate of WIMP and Ca is the specific annihilation rate related 
to the Sun volume. The second term in equation Q4.1|) has the meaning of the rate of 
WIMPs pair annihilations: 2 • Ta = CaN 2 . The process of WIMPs evaporation does 
not included in ( |4.1| ) since it is negligible for WIMP mass larger a few Ge V / c 2 [|30| , |32 



It is easy to show that when equilibrium time t eq = (Cc • Ca) is much less than 
the live time of the solar system the annihilation rate is the half of capture rate 

r A = -y -tanh^tsunAeq). (4.2) 

since the number of WIMPs is solved as 

N= \ 7T -tanh^sunAeq)- (4.3) 

V 



The quantity Ca is determined by following expression |30. p2| 



r _ (°"ann ' v) , ™ x ',3/2 (a „\ 

where (a ann ■ v) is the neutralino total annihilation cross section thermally averaged 
with the relative neutralino velocity v and Vs un = 5.7 ■ 10 27 cm 3 is the effective volume 
of the Sun. 

For calculation of the capture rate we use simplified model [[33] of the solar potential 



which can be parametrized by escape velocity v es (r) on distance r from the center of 
the Sun as follows 

vUr)=vl-^{vl-vl), (4.5) 
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where v c = 1354 km/s, v s = 795 km/s and M(r) is the mass of the Sun inside a sphere 
of radius r. In this approximation the capture rate can be written in an analytical 
form [0, H, ^ 



C 



M Q p x v 



{ (l + a)V 2 ' 



exp(— af] 2 ) 



{l + af/ 2 {A 2 c -A 2 s ) 



x (4.6) 



A, A 



+ 



1 1 + a 

2 ~ 
exp(— bfj 2 



, . . erf - erf (A) )+ — [ A_e" A + - A + e 
a — o / V / a/7t 



-.4- 



(a-6)(l + 6)V2(^2_ A 2 ) 



-(a-fe)A 2 



[2erf(f/) - erf + erf (A. 



-4c 



where rrii is the mass of z-th chemical element, /j is its mass content in the Sun, cij 
is total (the sum of SI and SD contributions) neutralino-nuclei elastic cross section, 



3^4 0) M 
2v 2 ix. 



A s = A\ v =v and A c = A\ Ves=Vc and the following notations are 



A(r) 
introduced 

A ± = A±fj, A± = A±fj, A = A{i + a) 1 / 2 , A = A{i + b) 1/2 , 



where En 



n 



3h 2 



n 



'l + aY/ 2 ' 



and R 



n 



v 



m x v 2 

/i + 1 



(4.7) 



0-91 (civ) 



1/3 



x 10 13 cm. Here we assume the 



2m N R 2 alLKX ±L ~ L"' i7i VGeW +0.3 

isothermal spherical model for halo WIMPs distribution, where r\ is equal to 1. Ex- 
pression ( [4.6[ ) is obtained with additional assumption the all chemical elements are 
distributed uniformly in the Sun. This allows to perform analytic integration over 
inner part of the Sun and results in the lengthy term in outer curly brackets. In 
our calculations we use uniformly distributed solar elements obtained from GS98 solar 
model presented in [^7 . 



Eq. ( f4.6|) shows explicitly the direct functional dependence of the capture rate on 
using approximations and values of scattering cross sections, velocity distribution, lo- 
cal dark matter density and solar model. In discussion of physical implications of the 
calculated capture rate (and, therefore, annihilation rate) one should bear in mind all 



related uncertainties. In particularly, according to recent analysis in Ref. ||35|| , error 



in determination of Ag in calculation of SD cross section results in 10% uncertanty 
in annihilation rate, while lack of exact knowledge of a n N can give rather substantial 
errors (up to order of magnitude) in annihilation rate if capture rate is saturated by 
SI interactions. The accuracy of using the approximation (|4.5| ) and assumption about 
uniformly distributed elements in the Sun with respect to e.g. AGSS09 solar pro- 
file |27j is within 11%. Also as it follows from cosmological N-body simulations [|36| the 
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Table 1: Benchmark NMSSM points. 



local dark matter density p l ° c might be about 20% larger than usually adopted value 
0.3 GeV/cm 3 . 

Now we consider the NMSSM benchmark points proposed and discussed in the 
paper [22| in the search strategies at the LHC In the Table [I] we present the physical 
quantities relevant for the dark matter searches in these models: the annihilation rates, 
the main annihilation channel and also SI and SD elastic cross sections. We note in 
passing that the models P1-P3 and P5 correspond to almost pure bino dark matter, 
while in the models P5 the lightest neutralino is singlino. Comparison with the latest 
results of the direct searches (see Sec. |5|) shows that even with account of the error bars 
(see next section for an illustrative figure) the points P4 and P5 are already excluded, 
while models P1-P3 can not be probed in the nearest future direct searches and at 
neutrino experiments. 

Annihilation rates V a calculated in described frames of Set 1 and Set 2 models are 
plotted in Fig.|l] as a function of neutralino mass. The upper limits Y u A ppLim obtained 
at 90% c.l. in neutrino experiments of the Baksan |J7| and the AMANDA [[38| collab- 
orations are also shown. In Fig.0 we markup points of the NMSSM models by colours 
correspondent to different dominant annihilation channels, with more than 0.5 branch- 
ing ratio in each case. There is obvious difference between these two SUSY samples 
from point of view annihilation branches. Annihilations to bb in mSUGRA-like model 
overcomes others for almost all neutralino masses. Consequently, these 66-like models 
can be excluded only by low energy threshold experiments as the Baksan |37[ |39], |4(| 
and the Super Kamiokande j|T| telescopes. A role of r + r~ channel is not remarkable 
in mSUGRA-like model for low neutralino masses, while we find expected competition 
of branching ratios of bb and t + t~ for Set 2. In mass range near 60 — 80 GeV/c 2 there 
is a decrease in the number of points due to closeness of possible value of the scalar 
Higgs mass to double neutralino mass. However the models inside this mass interval 
have a high potential to be seen or excluded by neutrino telescopes. The same conclu- 
sion can be drawn for scalars channels (i.e., annihilation of neutralinos to scalar and 
pseudoscalar Higgs bosons) for Set 2. 

In the frames of considered models the lightest neutralino could be singlino with 
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light mass even less than 50 GeV/c 2 . Not many singlinos satisfy condition Q x < 0.3 
and only small amount of them survive imposing 3cr WMAP density cut as it can be 
seen in Fig. |[ However singlinos appear to be among the most promising candidates 
in Set 2 in searching for neutralino annihilation signature in neutrino telescopes. To 
see which annihilation channel plays dominant role we show the rest of the models in 
Fig. ^ separately and respectively to colours in Fig. [TJ. 

We note that for neutralino mass less than 80 GeV/c 2 main annihilation channels 
are quark and lepton pairs. Neutralino annihilating mainly to scalars has larger anni- 
hilation rate then one with leptonic channel. As we see the Baksan exclusion line of 
YUppLim g oeg c i ose iy to model candidates. The AMANDA lower (solid) line of limit 
could be compared here only with r + r _ , W + W~ and tt points while for bb channels 
(green points) their Y^ ppLm limits [[38] (dashed line) are in two orders higher. 



5. Limits on neutralino- nucleon cross sections 



The basic conclusion followed from the equation (|4.1| ) is that annihilation rate is 
determined by WIMPs scattering cross sections off solar matter and in equilibrium does 
not depend on their annihilation cross sections 

r* = % (5.D 

In this case annihilation rate ( |4.2| ) can be divided into two pieces which correspond to 
different type of contribution (SI or SD) to neutralino-nuclei interactions 

T A = T A {a SI ) + r A {a SD ), (5.2) 



where Ta(&si) and Ta(ctsd) are determined by (|5.1j) and ( |4.6|) , where one should keep 
only SI or SD part of neutralino-nuclei cross section respectively. To demonstrate the 
relative importance of SI and SD contributions we plot their ratio in Fig. |] for Set 1 and 
Set 2 in dependence on neutralino mass and dominant annihilation channels (here and 
below we present in figures only the models within 3cr WMAP region for dark matter 
abundance). As expected the axial interactions play dominant role for relatively light 
neutralinos (m x < 200 — 500 GeV) while annihilation rate for heavier neutralino is 
saturated in main part by scalar interactions. Hence for large neutralino masses the 
uncertainties in the calculations of the annihilation rate increase. 



As it has been shown by G. Wikstrom and J. Edsjo [P3l , if the WIMP's processes 
in the Sun reach exact equilibrium one can derive upper limits on WIMP scattering 
cross sections from upper limits on upward going muon fluxes obtained by neutrino 
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telescopes. Consequently these results simplify comparison with direct WIMP searches 
via nuclear recoil reactions. 

Early in the Baksan neutralino search it was mentioned a preference of limits 
on annihilation rate Ta for the purpose of comparison with model prediction because 
muon flux $^ calculated for some model depends on details of the experiment, e.g. on 
its energy threshold. The magnitude of annihilation rate gives the scale for neutrino 
flux from a distant source (like the Sun at the distance R), that is 

Z^ B i-TET> = (5-3) 



dE v AnR 2 ^ dE, 



where dN„./dE v is the differential spectrum of j-th neutrinos at the surface of the 
the Sun produced in i-th annihilation channel. The branching ratios Bi depend on 
SUSY parameters as we show above and obviously drive the shape of the final neutrino 
spectrum, which could be either mainly "soft"-like from annihilation channels into two 
quarks (e.g. a pair of bb) or mainly "hard" -like from lepton decay channels of two 
bosons (e.g. a pair of iy + VU~). 

Definition of an upper limit on neutralino annihilation rate is read from the equa- 



tion 40 



tt t- ]\[U PP Lim,tfl\ i 

rU/P Um = W x 1 (5 4) 



aim 



where N UppLtm ( r d) is an experimental limit at a given confidence level on number of 
events and for a given half-angle of open cone i? toward the Sun, T is a time of observa- 
tion, and s$ is a fraction of total number of neutrino events collected within the angle i? 
from expected neutrino source. A value P an n is neutrino telescope probability to detect 
muon per one pair of neutralino annihilation. Here the upward going muon is produced 
in charged current neutrino interactions in the Earth. To get the probability P an n the 
detailed Monte Carlo simulations are required for the detector response with generator 
of neutrino energy spectrum ( |5.3j ) and further propagation of oscillating and interacting 
neutrinos in matter (see e.g. Refs. p9"| ) and also induced muons on the route to the 
detector. From these studies a relation between neutralino mass and the angle $ could 
be found. The latest results with 90% c.l. upper limits of annihilation rates in the Sun 



are presented by both the AMANDA and the IceCube |42| collaborations and as 
well by the Baksan Underground Scintillator Telescope in Ref. 

We suggest to implement obtained upper limits on annihilation rate Y^ ppLm to 
put upper limits on neutralino SD and SI scattering cross sections off protons (see also 
Eq.(15) in Ref. |23|) 

4^ m im x ) = X SD (m x ) • TT Um {m x ) (5.5) 
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Table 2: Coefficients X s D , X s 1 and the upper limits <Jg^^ ,m , (T ^ p ^ jim ' at 90% c.l. on elastic 
%-proton cross sections recalculated from the Baksan 90% c.l. upper limits on annihilation 
rates. 



Coefficients for these calculations are following 



U ppLim , 
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SD,p 
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r x A (as D y 



a 



SI,p 



r x A (*siY 



(5.6) 



(5.7) 



They can be obtained from formulae above with the assumption of dominance either 
spin-dependent or spin-independent interactions. Clearly, these coefficients in contrast 
with ratios ^/T^, cr SI>p /^^ and &sd,p/&[i presented in paper J23j, are universal in the 
sense that they do not depend on the experimental setup (e.g., on energy threshold) 
and on the model of the WIMP dark matter, although it still contains uncertainties of 
astrophysical values (see Eq. ([4. 6|) ) . Coefficients X SD and X SI as functions of neutralino 
mass m x are plotted in Fig.[5[ 

In Table |2] we present the upper limits on neutralino-proton elastic cross sections 
CTsi)p irn an d °'si P p im wn i cn we derive from the Baksan 90% c.l. upper limits on annihi- 
lation rate by use coefficients ( |5.7|) which are also shown in Table ^| The neutralino 



mass values are the same as in the Baksan final table 40 



Graphical view of the obtained limits on cross sections is shown in Fig. || Here for 
comparison (and as well to check our procedure) we show both the upper limits on SI 



and SD cross sections presented by the AMANDA experiment in [33] and our current 
recalculations from the AMANDA upper limits on annihilation rates ||38|| . We observe 
nice coincidence. Therefore coefficients X SD and X SI presented in Table |2] and in Fig.|| 
could be used in recalculation of annihilation rate limits for other neutrino telescopes 
(such as the Baikal NT-200 [46| and ANTARES [[|7j]) in wide range of neutralino masses. 

In Figs. |7| and || neutralino-proton SI and SD cross sections calculated for NMSSM 
models are shown with the strongest experimental upper limits o s Y^ im and (y V s ^^ tm 
obtained at 90% c.l. in direct detections (CDMS gO], XENON10 |43|', ZEPLIN-III |4 
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and PICASSO @) and neutrino telescopes (AMANDA ||], Baksan g7[ and Ice- 
Cube [|2[]). We implement the DMTools |48| database to plot listed upper limits. 
Although main part of the DAMA || allowed ranges are excluded by different direct 
experiments, these 3a regions with and without ion channeling are also shown as shad- 
owed. One can see that the SI limits from the neutrino detectors are able to close 
the DAMA results for neutralino mass around 80 — 100 GeV/c 2 . The Baksan upper 
limits on SD cross sections close the values covered by the DAMA for neutralino masses 
down to 12 GeV/c 2 . At the lightest mass the Baksan cr^^ %ra level is the same as the 



latest results of the ZEPLIN-III |4| and PICASSO [45J. Finally the Baksan limits 



on SD cross sections resolve the "gap" |23j between limits of the direct detections and 
predicted SUSY models, here the NMSSM models. 

In Figs. [7] and |8| we markup the different neutralino type by colours. Namely we 
identify the lightest neutralino as bino if Zq X > 0.5, as wino if Zq 2 > 0.5, as higgsino 
if Z$ 3 + Z$ 4 > 0.5 and as singlino if Z$ 5 > 0.5. For mSUGRA-type models (Set 1) 
the lightest neutralino is bino for almost entire parameter space. However for the 
Set 2 there are quite a lot of models in which dark matter can be wino, higgsino or 
singlino. Note that in Fig. [8] we drop all models which are excluded by the limits on 
SI neutralino-proton cross section shown in Figs. [7| with account of g^n uncertainty. 
Namely we plot only those points which are not excluded at least for some values of a^N 
from the interval 55 — 73 MeV. Also one can see that population of bino-like neutralinos 
in the NMSSM has better prospects to be found or excluded in experimental search by 
both direct methods and neutrino detections. Two benchmark points P4 and P5 are 
shown with error bars for SI cross sections. From SI limits in Figs. [7] one can conclude 
about exclusion of P4 and P5 models by direct detection experiments. Also the point 
P4 in Figs. is close to being excluded by SD limits from the neutrino telescopes. 



6. Summary and conclusions 

We discussed here the prospects of neutrino telescope searches for neutralino dark 
matter in the framework of the NMSSM. In our study we analyse two types of the 
NMSSM scenarios, namely, mSUGRA-like models and general NMSSM without uni- 
fication of the soft terms and perform a scan over the NMSSM parameter space with 
using NMSSMTools package. We have found that the dominant annihilation channels 
in the mSUGRA-like models are very different from those of the general NMSSM. In 
case of the mSUGRA models the lightest neutralino is mainly bino and it annihilates 
dominantly into bb pairs. However, neutralino annihilating into W + W~ and ti pairs 
is also very promising for searches at neutrino telescopes. For general NMSSM model 
we have found expected competition between bb and r + r~ channels of annihilation for 
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neutralino mass less than mass of W-boson, while for heavier neutralino W + W~ and tt 
channels are comparable with those into scalar pairs. We found that NMSSM neutrali- 
nos with mass less about 200 GeV/c 2 scatter off nucleons in the Sun mainly by axial 
interactions and with increasing of the mass value the part of the scalar interactions also 
increases such that at masses higher than 500 GeV/c 2 it is almost spin-independent 
scattering. In general, for both types of NMSSM models singlino and bino are the 
most probable models to be tested either at the neutrino telescopes with low energy 
threshold through neutralino axial-vector (SD) interactions in the Sun or in the direct 
detection experiments sensible to scalar (SI) neutralino interactions. 

Depending on neutralino properties there are essentially different prospects in 
search for them at neutrino telescopes with low and high energy thresholds. The 
type of annihilation products determines fractions of low (soft) and high (hard) energy 
spectra of outcoming particles including neutrinos. Probability of neutrino telescope 
to detect neutrinos from predicted neutralino models absorbs dependence on energy 
threshold of the detector for neutrino induced upward going muons. Presently among 
the operating neutrino telescopes the Baksan Underground Scintillator Telescope and 
the Super Kamiokande have low energy thresholds (about 1 GeV). The Baksan 90% 
c.l. upper limits on annihilation rate in the Sun result from absence neutralino signal in 
all expected mass range from dozen GeV/c 2 . The large detectors with high thresholds 
e.g. AMANDA or IceCube present the upper limits on annihilation rates for separated 
annihilation branches with soft and hard neutrino spectra. Therefore all limits can be 
compared only in context of a particular neutralino annihilation channel. 

Finally, we have discussed a relation between the limits on neutralino annihilation 
rate in the Sun and the limits on neutralino-proton SI and SD cross section. Under 
assumption of equilibrium of the neutralino processes in the Sun we derive 90% c.l. 
upper limits on SI and SD neutralino-proton cross sections from the upper limits on 
annihilation rates for the Baksan Underground Scintillator Telescope. Also we tabulate 
the relevant coefficients for calculation the above limits, since they are independent on 
the details of a particular experiment and the type of WIMP dark matter model. 
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Figure 1: Annihilation rates vs neutralino mass for Set 1 (up) and Set 2 (bottom) of 
parameters. Here all points are in the interval < Vt x < 0.3. We show the dominant 
annihilation channel (Br > 0.5) by different colours. 
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Figure 2: Annihilation rates vs neutralino mass for Set 1 (up) and Set 2 (bottom) of 
parameters. 
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Figure 3: Annihilation rates vs neutralino mass for Set 1 (up) and Set 2 (bottom) of 
parameters. Here all points are in 3a WMAP interval: 0.0913 < Qdm < 0.1285. We show 
the dominant annihilation channel (Br > 0.5) by different colours. 
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Figure 4: Comparison of relative contribution of spin-dependent and spin-independent 
interactions to annihilation rate. All conventions are the same as in Fig. 0. 
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Figure 5: Coefficients Xsi(m x ) (up) and Xsd(^x) (bottom). 
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Figure 6: The upper limits at 90% c.l. on neutralino-proton spin-independent (up) and 
spin-dependent (bottom) cross sections for the Baksan underground scintillator telescope 
recalculated from their 90% c.l. limits on annihilation rates (37]] (red lines). Also shown 
are the AMANDA upper limits at 90% c.l. on SI and SD cross sections presented in 3£] 
(triangles) and our recalculated values from the AMANDA upper limits on annihilation rates 
(blue lines). 



- 25 - 




Neutralino mass [GeV/c 2 ] 




Neutralino mass [GeV/c 2 ] 



Figure 7: Neutralino-proton spin-independent cross section vs neutralino mass for Set 1 
(up) and Set 2 (bottom) of parameters and experimental limits at 90% c.l. Neutralino content 
is shown by different colours. 
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Figure 8: Neutralino-proton spin-dependent cross section vs neutralino mass for Set 1 (up) 
and Set 2 (bottom) of parameters and experimental limits. The notations are as in Fig.0. All 
points are constraint by SI cross section experimental limits at 90% c.l. 
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